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We propose that in a system with pocket Fermi surfaces, a pairing state with a finite total 
momentum q tot like the Fulde-Ferrell-Larkin-Ovchinnikov state can be stabilized even without 
a magnetic field. When a pair is composed of electrons on a pocket Fermi surface whose center 
is not located at F point, the pair inevitably has finite q tot . To investigate this possibility, we 
consider a two-orbital model on a square lattice that can realize pocket Fermi surfaces and we 
apply fluctuation exchange approximation. Then, by changing the electron number n per site, 
we indeed find that such superconducting states with finite q tot are stabilized when the system 
has pocket Fermi surfaces. 

KEYWORDS: FFLO state, pocket Fermi surfaces, orbital degree of freedom, fluctuation exchange approx- 
imation 



The pairing states of superconductivity are classified 
on the basis of symmetries that satisfy the Pauli prin- 
ciple. *) Among the pairing states, the s-wave state is 
the simplest and is called conventional superconductivity. 
The other superconducting states are anisotropic (uncon- 
ventional) ones. We can categorize the superconducting 
states into spin-singlet and spin-triplet states. The spin- 
singlet states have even parity, e.g., s-wave, under the 
inversion of space due to the Pauli principle. The spin- 
triplet states have odd parity, and in these states, the 
spin degree of freedom is active and many phases can 
occur in a system. 

If we can introduce an additional degree of freedom 
into a system, we can find exotic superconducting states 
that are beyond the ordinary categorization. For this pur- 
pose, the superconductivity in systems with an orbital 
degree of freedom, which is not considered in the above 
discussion on superconducting symmetry, has been stud- 
ied recently. For example, the effects of the orbital degree 
of freedom on superconductivity have been studied for a 
two-orbital Hubbard model with the same dispersion for 
both orbitals, by a mean-field theory, 2 > dynamical mean- 
field theory, 3,4 ) and a fluctuation exchange (FLEX) ap- 
proximation. 5 ) These studies have revealed that an s- 
wave spin-triplet state and a p-wave spin-singlet state, 
which satisfy the Pauli principle by composing the or- 
bital state of a pair antisymmetrically, can be stabilized 
in the two-orbital Hubbard model. 

Although the two-orbital Hubbard model provides 
such an interesting possibility of superconductivity in a 
multi-orbital system, it describes a system with a rather 
high orbital symmetry, that is, the Fermi surfaces of these 
orbitals are exactly the same. To discuss more realistic 
situations, we should improve the two-orbital Hubbard 
model. In particular, we should include the effects of or- 
bital anisotropy, for example, orbital-dependent hopping 
integrals and the transformation property of orbitals. 
These properties are important for magnetism in d- 
electron systems 6, 7 ' and /-electron systems, 7 ~ 9 > and they 
should also be important for superconductivity. Among 
the above-mentioned properties, the multi-Fermi-surface 



Fig. 1. Pocket Fermi surfaces. A pair with a total momentum q tot 
is composed of electrons with fc and — fc + q tot which are on the 
same Fermi surface. This pair can be scattered to another Fermi 
surface. 



nature due to the orbital-dependent hopping is impor- 
tant, since a pair on a Fermi surface can interact with a 
pair on another Fermi surface, which may provide a new 
pairing mechanism. 

Among possible multi-Fermi-surface cases, a system 
with pocket Fermi surfaces, as shown in Fig. 1, is in- 
teresting for superconductivity. In such a system, a pair 
on a Fermi surface has a finite total momentum q tot like 
the Fuldc-Fcrrell-Larkin-Ovchinnikov (FFLO) state 10, ir > 
even without a magnetic field. When the electrons com- 
prised in this pair are scattered to another Fermi surface 
by fluctuations such as spin fluctuations and can form 
a pair on that Fermi surface, such a pairing state with 
finite q tot can be stabilized. 

In this Letter, in order to explicitly discuss the pos- 
sibility of such exotic superconductivity, we consider an 
e g orbital model on a square lattice as an example. This 
model can easily realize pocket Fermi surfaces even if 
we consider only nearest-neighbor hopping. To investi- 
gate the possible superconducting states in this model, 
we apply FLEX approximation that has been extended 
to multi-orbital models. 5, 12 ~ 16 ) Then, we show that such 
pairing states with a finite total momentum like the 
FFLO state are stabilized in a system with pocket Fermi 
surfaces. 



J. Phys. Soc. Jpn. 



Letter 



Katsunori KUBO 



71 

-71 



n=0.5 







n=0.75 



n=l 




n=1.25 



«=1.5 



n=1.75 




«=1.9 



-7C 7C-7C 7C 



J* 



71 



*0 



-7C 



— i 1 1 — 

o o 

O Q 

I I I 



-7C 7C 



Fig. 2. Fermi surfaces of the model for (dda) = — (ddS) for several 
electron numbers n per site. 



We consider a tight-binding model for e g orbitals given 



by 



H 






£kTT>cl T(T c kT >a + U'^njT^njri 



where Ci Ta is the annihilation operator of the electron 
at site i with orbital r (= 1 or 2) and spin a (=j or 
I), Ckra is its Fourier transform, rii Ta = cj T0 .Ci T(T , and 
n iT = ^2 ni T cr- The coupling constants U, V, J, and J' 
denote the intra-orbital Coulomb, inter-orbital Coulomb, 
exchange, and pair- hopping interactions, respectively. In 
this study, we use the relations U = U' + J + J' and 
J = J'. 17 ) 

We consider the hopping integrals for the nearest- 
neighbor sites on a square lattice, and the coefficients 
of the kinetic energy terms in Eq. (1) are generally given 
by Cfen = \[d(ddcr) + (dd5)](cosk x + cosfc y ), e fc22 = 
\[{dda) + 3(dd5)](cos k x + cosk y ), and e fc i 2 = e k2 \ = 

— ^[(dda) — (ddS)](cosk x ~ cosk y ), where (dda) and 
(ddS) are Slater-Koster integrals, and the lattice constant 
is set to unity. 18 ) 

To realize pocket Fermi surfaces in a system such as 
that shown in Fig. 1, we set (dda) = —(ddS) = —?=t. 
Figure 2 shows the Fermi surfaces for several electron 



numbers n per site. Due to the electron-hole symmetry, 
it is enough to consider < n < 2. For n — 2, the Fermi 
surfaces disappear in this model. Pocket Fermi surfaces 
are realized for n > 1.25. In this system, if an electron 
with — k is on a Fermi surface, the electron with —k + Q 
is on another Fermi surface where Q = (tt, tt), as shown 
in Fig. 2 for n = 1.75 as an example. Then, it is possi- 
ble to form a pair with the total momentum q tot — Q by 
electrons with k and — k + Q. Thus, in this study, we con- 
sider superconducting states with q tot = Q in addition 
to the ordinary superconducting states with q tot = (0, 0). 
The Green's function for the present two-orbital model 
is expressed by a 2 x 2 matrix. In the normal phase, 
the Dyson-Gorkov equation for the Green's function in 
matrix form is given by 



G{k) = G (0) (fc) + G (0) (fc)S(fc)G(fc), 



(2) 



where k = (fe, ie„) and e n = (2n + 1)ttT is the Matsubara 
frequency for fermions with an integer n and a tempera- 
ture T. The non-interacting Green's function is given by 
G(°)(fc,ie„) = [ie„ — e k + /i] _1 , where /x is the chemical 
potential. In the FLEX approximation, the self-energy is 
given by 



T 



^rr>{k) = - ]T V, 



TT\ ;t't2 



(q)G TlT2 (k - q), (3) 



Q,T~1,T2 



where TV is the number of lattice sites, q = (q, iw m ), and 
u> m = 2rrnrT is the Matsubara frequency for bosons with 
an integer m. The matrix V(q) is written as 



V(q) =^[U s X s (q)U s - U s X ^(q)U s /2 + U s ] 
+ \[U c X c (<l)U c - U c x m {q)U c /2 - U% 



(4) 



where the matrix elements of U s and U c are given by 

TTS _ TJS _ TTC _ TTC _ TT TTS _ 

u ll;ll — u 22;22 — ^ll;!! " ^22:22 ~ u > Ly ll;22 — 
TTS _ j ttc _ tjc _ 977/ _ J TTS _ 77s _ 

t7 22;ll — J i c 'll;22 _ u 22:ll _ Lu J 1 L/ 12:12 _ c; 21;21 _ 

U , C/l2;12 = ^21;21 = — U + 2J, C'l2;21 = ^21;12 = 

^12-21 = ^2112 = "^'i the other matrix elements are 
zero. The susceptibilities X s (g) and X c (<z) are given by 
X s (9) = X (( %)[1 - ^l?)]' 1 and x c (<7) = X (0) (<z)[l + 
U c \ (q)]^ 1 j respectively, in matrix forms. The ma- 
trix elements of x (?) are defined by Xtit 2 :t3t 4 (<z) = 
-|f Efe G^rs (k + q)G TlT2 (k). We solve Eqs. (2)-(4) self- 
consistcntly. Then, we can calculate response functions; 
for example, the spin susceptibility in the FLEX approx- 
imation is given by x(q) = ^Y, t .t' X s tt;tV (9)- 

The linearized gap equation for the anomalous self- 
energy 4> TT '{k; q tot ) is given by 



T 



0rr'( fc ;9tot) = -T7 J2 



V r \ i]T2Tl (k-k>) 



fe / ) Tl,T2,T3,T4 



(5) 



xG TlT3 (k')(j) T (k';q tot )G T2Ti (-k' + g tat , -ie„), 



where the spin state is denoted by ^ — singlet or triplet, 
and q tot denotes the total momentum of the pair. The 
superconducting transition temperature is given by the 
temperature for which Eq. (5) has a nontrivial solution. 
The effective pairing interactions V^ (q) are written as 



F singlet (g) = [w s {q) _ V c {q)]/2 , 



(6) 
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V t ^ let (q) = [-V s (q)-V c (q)}/2, 



(7) 



where V s (q) = [U s x s (q)U s + U s /2] and V c (q) = 

[u c x c (q)u c -u c /2]. 

Before presenting the calculation results, we discuss 
plausible candidates for the pairing symmetry from 
Eqs. (5)— (7). First, in this discussion, we ignore the or- 
bital indices in the above equations for simplicity. In this 
model, we find that only the spin susceptibility becomes 
large. Then, we expect an anisotropic superconducting 
state for the spin-singlet channel, since V slnglet (q) be- 
comes large and has a positive value and the anoma- 
lous self-energy has to change its sign in the k space, 
as understood from Eq. (5). For the spin-triplet channel, 
an s-wave state is favorable, since V tTlplet (q) becomes 
negative and the anomalous self-energy does not have 
to change its sign. The above consideration is appropri- 
ate for orbital-parallel components, i.e., (^(fc; <7 t ot) an d 
4>22{k', 9tot)i as l° n g as the Pauli principle is satisfied. 
On the other hand, for orbital-antiparallel components 
(j)\ 2 (k; Qtot) an d 0li (^j 9tot)j orbital indices transform fol- 
lowing the d x 2_ y 2 symmetry under symmetry operations 
of the lattice, and thus, the anomalous self-energy be- 
longing to the d x 2_ y 2 symmetry does not have to change 
its sign in the k space. In other words, the symmetry 
in the k space of a d x 2_ y 2-wa,ve state for the orbital- 
antiparallel components is the same as that of an s-wave 
state for the orbital-parallel components. Thus, spin- 
triplet d a ,2_ J/ 2-wave states and spin-singlet states other 
than those with the d x 2_ y 2 symmetry are favorable for 
the orbital-antiparallel components. Here, we note that 
orbital states mix in general. However, it is expected 
that when the Hund's rule coupling J is large and the 
inter-orbital Coulomb interaction U' is small, the orbital- 
antiparallel components become dominant. 

Here, we show the results for a 64 x 64 lattice. In this 
study, we fix the value of the intra-orbital Coulomb inter- 
action U = 6t and vary J (= J'). Then, the inter-orbital 
Coulomb interaction is given by U' = U—2J. In Fig. 3(a), 
we show the n dependence of the static spin susceptibil- 
ity X(q) = X(<?,iw m = 0) at q = q max for J = 0, t, and 
2i at T — 0.005£, where g max is defined as the wave vec- 
tor at which x(l) becomes the largest. For comparison 
purposes, we also show x(<Z max ) for the non-interacting 
system. The spin susceptibility is enhanced by the inter- 
actions, in particular, by the Hund's rule coupling. On 
the other hand, q max does not depend so much on the in- 
teractions as shown in Fig. 3(b). This fact indicates that 
the characteristic wave vector q max is almost determined 
by the property of the non- interacting system, i.e., the 
Fermi-surface structure. 

As shown in Fig. 3(b), the characteristic wave vec- 
tor q max takes various values depending on n around 
n = 0.7-1.5, since the Fermi surface structure changes 
around this area, as shown in Fig. 2, e.g, the Fermi sur- 
faces consist of electron sheets for n = 1, while they 
consist of hole pockets for n = 1.25. Thus, several fluc- 
tuations exist in this region, and a d x 2_ y 2--wave spin- 
triplet state with dominant orbital-antiparallel compo- 
nents may appear if such fluctuations in a wide q region 
are available. However, around n — 0.7—1.1, the spin fluc- 
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Fig. 3. (a) spin susceptibility x(q) at q = q max and (b) qr max for 
J = (circles), t (triangles), and 2t (squares) as functions of n 
for U = 6t at T = 0.005i. The thin lines are for U = 0. 
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Fig. 4. Fermi surfaces for n = 1.75. 



tuations are very strong for J = 2t and thus magnetic 
order should occur. For n > 1.5, we expect the appear- 
ance of a p-wave spin-singlet state because of the follow- 
ing reason. For the purpose of explanation, we show the 
Fermi surfaces for n = 1.75 in Fig. 4. The characteristic 
wave vector q max connects electrons at k = q max /2 and 
—k. Thus, electrons around k and k + Q comprising a 
pair on a Fermi surface can be scattered to another Fermi 
surface by the spin fluctuations and can comprise a pair 
with momenta around — k and k + Q. These pairs are the 
space inversions of each other, and to utilize such fluc- 
tuations, the p-wave spin-singlet state is favorable since 
the anomalous self-energy changes its sign under the in- 
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Fig. 5. (a) Transition temperatures for U = 6t. The solid trian- 
gles and squares denote antiferromagnetic (AF) transition tem- 
peratures for J = t and 2i, respectively. The open circles and 
triangles denote superconducting transition temperatures for the 
d x 2_ y 2-wave spin-triplet state and p-wave spin-singlet state, re- 
spectively, for J = 2t. In both superconducting states, Cooper 
pairs have the total momentum Q = (tt,tt). (b) Transition tem- 
peratures around n = 1.6 for J = 2t. 



version k — ► —k. These d x 2_ y 2-wave spin-triplet and p- 
wave spin-singlet states should be dominated by orbital- 
antiparallcl components because of the Pauli principle. 
The orbital-parallel components should be odd in fre- 
quency for these superconducting states and play only a 
minor role. 

Figure 5(a) shows the highest transition temperatures 
among transition temperatures for all superconducting 
states allowed in the model and for antiferromagnetic 
states as functions of n for J = t and 2i. Figure 5(b) 
shows the highest transition temperatures for J = It 
around n — 1.6. The antiferromagnetic transition tem- 
peratures are determined by x(9max) — 70/t, and the 
system is always in the normal phase for J = at 
T > 0.005i. The antiferromagnetic phase extends by in- 
creasing the Hund's rule coupling. With regard to super- 
conductivity, we cannot find any superconducting phase 
within T > 0.005i for J = t. For J = 2i, a p-wave spin- 
singlet state with q tot — Q appears, as is expected from 



Fig. 4, at n ~ 1.75 where the antiferromagnetic transition 
temperature tends to zero. The d x 2_ y 2-wave spin-triplet 
state with q tot = Q appears in a wider region n ~ 1.3- 
1.7, since the dj;2_ y 2-wave spin-triplet state is insensi- 
tive to the value of the characteristic wave vector q max , 
similar to the s-wave state in the two-orbital Hubbard 
model. 5 ) It should be noted that these superconducting 
states with finite q tot appear only in the region where 
the system has pocket Fermi surfaces. 

To summarize, we have pointed out that a pairing state 
with a finite total momentum like the FFLO state is pos- 
sible for a system with pocket Fermi surfaces. As an ex- 
ample, we have studied a model for e g orbitals by apply- 
ing FLEX approximation. Then, we have found that the 
p-wave spin-singlet and d x 2_ y 2-wave spin-triplet states 
with q tot = Q = (-7T, 7r) are realized at electron number n 
where the system has pocket Fermi surfaces. The mech- 
anism to stabilize the superconducting states with finite 
q tot can be applied to systems other than the e g orbital 
system as long as such pocket Fermi surfaces exist. 
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